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Abstract 

We generalize an already proposed protocol for quantum state transfer to 
spin chains of arbitrary spin. An arbitrary unknown d— level state is trans- 
ferred through a chain with rather good fidelity by the natural dynamics of the 
chain. We compare the performance of this protocol for various values of d. 
A by-product of our study is a much simpler method for picking up the state 
at the destination as compared with the one proposed previously. We also dis- 
cuss entanglement distribution through such chains and show that the quality of 
entanglement transition increases with the number of levels d. 
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1 Introduction 



Since the proposal of S. Bose [lj for transferring quantum states via natural evolution 
of quantum spin one-half chains, there have been many types of extensions of this 
idea in various directions. For example it has been shown that perfect transfer is 
possible for a special class of Hamiltonians, called mirror-periodic [2JE]. It has also 
been shown that one can achieve better fidelities either by using multiple chains [I] 
or by allowing the parties to have access to more than one site of the chain [5] or 
by using chains with longer range interaction than nearest neighbor [6j. The effect 
of thermal fluctuations [7] and decoherence [5J [9] have also been taken into account. 
Some other aspects of this protocol have been studied in [TOl \TT[ [T2]. 
However to our knowledge there has been no attempt to generalize this proposal to 
chains of particles of arbitrary spin. The aim of this paper is to extend this proposal 
in this new and fundamental direction. There are good reasons why such an exten- 
sion is worthwhile. First, until a particular experimental proposal for qubit quantum 
computer is widely accepted as the platform for implementation of quantum com- 
puters, we have to formulate various theoretical protocols for particles with arbitrary 
number of levels, the so-called qudits. In fact for this reason, various protocols of 
quantum computation and information, like cloning [131 114j . cryptography [15] and 
teleportation [161 [T7] have been generalized to d- dimensional systems. Second, from 
purely theoretical point of view we will learn very much in developing a particular 
scheme like quantum state transfer in a way such that the role of dimensionality can 
be studied in detail. In fact the work of Bose PQ can be rephrased in a way which 
demands such an extension in a quite natural way: It is well known that a quantum 
state can be transferred perfectly through a chain by sequential application of the 
swap operator defined as P\a,j3) = \P,a). However this method requires control on 
every qubit throughout the chain. Instead in [1] a state is coupled to left hand side 
of a spin one-half chain, governed by a ferromagnetic Heisenberg Hamiltonian 



where Si = ^a{, cij's are the Pauli operators, J > is the coupling constant and B 
is the magnetic field. Then the natural evolution of this chain will transfer the state 
to the right hand side, with a good fidelity provided that the state is extracted at an 
optimal time. This method can be named random swapping of a state. The reason 
is that using the identity 



where P is the permutation operator (P\a,(3) = \/3,a)), and a suitable redefinition of 
constants, H can be rewritten as 
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On the sector with fixed total spin S z , the evolution operator is equivalent to U = 
g'^E.^.Hi^ where we have set fi = 1 and ignored an overall phase. Thus for an 
infinitesimal time step e, we have 

m + e )> = m)) + E ijep iti+ im)), 

i 

which shows that the state \i()(t+e)) is obtained by adding to \ip(t)} an equal superpo- 
sition of states in which the spins of two adjacent sites have been swapped, hence the 
name random swapping. Thus the result of [1] can be rephrased in the following form: 
for qubits, random swapping achieves a fidelity which is reasonably good compared to 
that of sequential swapping (for which U = ]T, Pa+i)- In particular when the length 
of the chain is 4, the results of pQ imply that sequential and random swapping attain 
almost equal fidelity. Once interpreted in this way, we can ask naturally what form 
this comparison takes for states of arbitrary dimensions. 

We should note that the Hamiltonian ([1]) can always be expressed in terms of nearest- 
neighbor scalar spin interaction terms although in each dimension it takes a specific 
form, for example in dimension d = 3, it takes the form 

JV 

H = — J E (Si • Sj+i + (Sj • Sj + i) 2 ) + B E Siz- 

i i=l 

We will find that for a fixed distance, the fidelity decreases with dimension d, but 
reaches a saturated value depending on the distance and when the sender and the 
receiver are 4 sites apart, nearly perfect transfer is possible for any dimension d. As 
a by-product of our study, we will propose a much simpler method for state transfer, 
one in which the magnetic field is kept to a vanishingly small value, instead of tuning 
it to a distance-dependent value as in the original protocol of pp. The structure of 
this paper is as follows. In section [2] we introduce the basic protocol in d dimensions, 
and derive the basic relations that we need in the sequel. In section [3] we study the 
problem of entanglement distribution in such chains. In section U] we conclude with a 
discussion. 



2 Quantum state transfer in chains of qudits 

Originally the problem of state transferring was considered for an open chain [1 J. How- 
ever in that same work it was shown that in a ring of size 2N one can as efficiently 
transfer states as in an open chain as long as the distance between the sender and the 
receiver is not longer than N. To use the the advantage of simplicity of eigenfunc- 
tions of the Hamiltonian, we consider a periodic chain of length N, where each site 
comprises a state of a d level system with basis states fj, = 0, 1, • • • , d — 1. The 
evolution of the chain is governed by the Hamiltonian, 

_ j N N 
# = ^-E( P M+i-l) + £E^> (2) 
i=l 4=1 
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where the operator P^i+i is the permutation operator on sites i and i + and S{ z is a 
diagonal operator acting on the states of site i as, S z \fi) = fj>\fj>), for fj, = 0,1, ...,d— 1. 
Note that S* 2 , when shifted suitably, plays the role of the third component of the spin 
operator. Thus B plays the role of a magnetic field in the ^direction. The Hamil- 
tonian reduces to the Heisenberg Hamiltonian for spin 1/2 states, and to the 
bilinear-biquadratic hamiltonian for spin 1. For other spins it contains high power of 
the term (Sj • Sj+i). We assume that B is positive. 

The ground state of this Hamiltonian is given by |0) = |0)® with energy E g = 0. 
The reason is the following. Since the permutation operator has the property P 2 = I, 
its eigenvalues are ±1, and the operator J(l — P^j+i) will be a positive operator with 
eigenvalues and 2J. Therefore in the absence of magnetic field, the Hamiltonian, 
being a sum of positive operators, is positive and since the states fJ, = 0, • ■ • d— 1 
all have zero energy, they form the degenerate ground state of H{B = 0). The mag- 
netic field only removes the degeneracy and lowers the energy of the state |0)® , with 
respect to others (note that in our notation |0) has the lowest value of spin compo- 
nent.) Since the Hamiltonian commutes with S z , and H{B = 0) can be diagonlaized 
in sectors with fixed z component of spin, this argument is valid for all values of the 
magnetic field B. 

We should stress that in the absence of magnetic field, the phase diagram (i.e. the 
character and long range order in the ground state) of ([2]), may be quite complicated. 
This will then affects crucially the quality of state transfer in such chains, a problem 
which has been recently investigated for spin 1 chains in |18j . In the presence of 
magnetic field however, the ground state has a simple ferromagnetic order given by 
the ground state \0)® N . 

Let us denote a state in which the i-th site has been exited to the level \x by 

i.e. 

\ta) = |o,---o,//,o,---o). 

The permutation operators in H only displace this state through the chain and hence 
the Hamiltonian can be diagonalized in each sector in which the number and type of 
excited local states is fixed. This is a consequence of a number of conservation laws, 
namely 

TV 

[H,Q m ]=0 , :=J2(Sz,i) m 

i=l 

for m = 1, 2, • • • , d — 1. In d = 2 dimensions only the charge is conserved. These 
conservation laws imply for example that a state like, 1 1, 1, - - - 0, 0, 0) can not evolve 
to a state like |2, 0, • • • , 0, 0, 0), since although their charge are equal they have 
different charges. 

The states with only one site excited are called one particle states and the subspace 
spanned by these vectors comprise the one-particle sector of the full Hilbert space. 
Let us denote by the one particle sector with charge equal to \i. The whole 
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one particle sector is 



^^Weyf)®...^- 1 ' 



The Hamiltonian can be diagonalized in vj^ with eigenvectors given by, 

1 N 

\K) = ~7W E eT l^), m = 1, 2, TV, 
ViV fc=i 

with energy given by = J — Jcos(^j^) + Bfi. For quantum state transferring we 
can consider site s as the sender of the system and site r as a receiver. The initial 
state that should be sent is 

d-l 

So the initial state of the system (the site s plus the chain) is, 

d-i 

|Y>(0)) = |V> S ) ® |0) = a |0) + E a> s >. 

In view of the fact that H\0) = 0, the state at time t will be, 

d-i N d-i 

m)) = «o|0) + E W~ iHt \Vs) = «o|0) + E E fksa^e- iB ^ k ). 

fi=l k=l /i=l 

In deriving this formula we have used the fact that [S Z ,H := = and the 

conservation laws which restricts the evolution to the one particle sector of fixed 
charges. We have also defined 

fks ■= (Vk\e~ iHt \Vs), 

which is indeed independent of /i and hence can be taken outside the sum. 

The state of site r which is acting as the receiver will be generally mixed, so is 
denoted by p r (t) and is obtained by tracing out the other sites. 

Pr (t) = tr?m))W)\ = (1 - E WU M s | 2 )|0><0| + E aoa*/£|0>M 

fi=i fi=i 

fi=l M,^=l 

Rearrangement of the right hand side yields 

Pr(t) = (l-P)\0){0\+P\<f>){<f>\, 

where 

P = |/ rs | 2 (l-|a | 2 ) + |a | 2 , 
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and 

y p n=l 

Alternatively we can say that the input state p s (0) is mapped to the output state 
p r (t) by the positive map 

M 

where the so-called Kraus operators are given by 



A = \0)(0\ + Y,f& MM, A^ = v /l-|W |0)(m|, M = l,-,rf-l- (3) 

n=i 

The fidelity between the received state p r (t) and the initial state p s (0) = \ip s )(' l Ps\ is 
defined by F = \{ip s \p r (t)\ip s )\ 2 which turns out to be 

d-l 

F(t) = |a | 2 + ]T|ao| 2 |cg 2 {i 

In the sequel we should maximize this fidelity when it is uniformly averaged over the 
input states. The average is defined by 

(F(t)) = J F(t)dU, 

where dU is an invariant (Haar) measure over the SU(d) group, normalized such that 
J dll = 1. The reason for this choice of measure is as follows. Let us fix a basis, 
like {|0), |1), • • • \d — 1)}. We take a fixed reference state like |0) and note that every 
arbitrary state can be obtained from |0) by the action of a unitary operator U, 
i.e. = C/|0), for some non-unique U G SU(d). However, any two unitary matrices 
U and Ug, where g £ SU(d — 1) leaves |0) invariant, lead to the same state 
Therefore a proper measure that prevents this multiple counting is a measure over 
U{d)/U{d — 1). However since every state is multiply counted equally ( by a factor 
which is exactly the volume of the group SU(d — 1)), this does not affect the final 
averaging and we can use the simple measure over U(d). Invariance of this measure 
under left multiplication, i.e. dU = d(gll) guarantees uniformity of the measure over 
the space of all states. In two dimensions one can avoid multiple counting in a simple 
way, since in this case U(2)/U(l) ~ SO(3)/SO(2) ~ S2 and therefore, one can use 
the measure over the 2 dimensional (Bloch) sphere to count every state once. This is 
the measure used by Bose in pp. 

For a d— dimensional normalized state = X^=o a ^\^) an invariant measure yields 
trivially (|a^| 2 ) = \ V/i. To calculate the other averages we use = C/|0) and write 

<M 4 > = (l«o| 4 ) = / Woo\ A dU = -^Yy V n, 



W£l 2 }+ E K\ 2 K\ 2 fMi 

u,v=\ 
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where for the last equality we have used a result from [19] on invariant integration 
on unitary groups. We can now calculate (|a^| 2 |a^| 2 ) for fj, ^ v. In view of the 
normalization of the state, we have 

|2|„ |2\_V^/|„ |4\ | |2i |2\ 



i = E<KI 2 M 2 > = EKI 4 > + E<K 

= + +d{d-i)(\a lt \ 2 \a v \ 2 ), 

thus we find ^ 

Using these results we can now calculate the average of fidelity over a uniform ensem- 
ble of input states. 

In order to write (F(t)) in a simple form we note that 

i E m = i E e-^ Bt f rs \ = \f rs \r d (Bt), 



where 



and 



(d-l)Bt 



sm ^ 

T d (Bt) 1 



Bt ' 



sin — 



£(/& + /*«) = 2cos( 7 r fl - ^)T d (Bt), 



where / rs = |/rs|e i7rs . Inserting these in ([5]) we find that 

For ci = 2 we recover the formula of pQ. 

It remains to calculate the explicit expression for the amplitudes f{? s . We note 
that 

N . -iJt N-l 

f rs = <Mr|e-^V.> = E(^l e "^W>(^l/*.> = V^ e^^e 1 ^ 1 . 

k=l k=0 

For large N, a closed formula for f rs can be obtained by writing the right hand side 
as an integral. Thus 

f rs ~ e — [** e iJtcos(8)+i8(r-s) = e -*(J*-|(r-.)) J r _ s (J t ) j Large jy, 
27T 7n 
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where J n is the Bessel function of the first kind of order n. It now remains to follow 
a definite strategy for picking up the state at the destination point r. As is clear 
from equation (0), there is a distinctive difference between dimension d = 2 and any 
other dimension, since in d = 2 we have T2(Bt) = 1 and the magnetic field enters 
only in one single term, namely the argument of cosine function. (Note that we can 
always re-scale the other coupling constant J to 1.) Thus in d = 2 there is rather a 
unique strategy, first suggested in pQ: at any given time t one finds the magnitude 
of B = Bit) which maximizes the cosine function to unity and then searches among 
the values of time t to determine the optimal time t op t for picking up the state. 
This will then determine the optimal value of the magnetic field through the relation 
B opt = B(t op t). Note that since t op t depends on the distance r — s between the sender 
and the receiver, the optimal value of the magnetic field also depends on this distance. 
This is an inconvenient feature of this strategy. 

A by-product of the present work is that a much simpler strategy can be used, 
namely: tune the magnetic field to a vanishingly small value, then the optimal time 
for picking up the state is almost independent of the magnetic field. To see this we 
note that for higher values of d, the magnetic field enters in two different ways in the 
final formula for the average fidelity, namely in the argument of the cosine function 
as in d = 2 and in the function T^(Bt). These two functions may have incompatible 
properties so that they may not be maximized simultaneously. The latter function 
is maximized when its argument Bt goes to zero, while the former function has a 
complicated dependence on t and B separately. Thus one can follow two different 
strategies for picking up the states at point r, the first one is exactly the same as in 
PQ, explained above. We can also follow a second much simpler strategy, which has 
the advantage of no need for distance-dependent tuning of the magnetic field. We 
simply apply a vanishingly small magnetic field (Bt « 1) for all the times involved 
in the transfer process. On the other hand B should not be vanishing so that we 
have a unique ferromagnetic ground state. This maximizes the function T^(Bt) to 
~ (d — 1). We are now left with a function which is entirely a function of t and can 
find for any distance r — s, the optimal time and the maximum average fidelity. 
In figure (1), we show the average fidelity for transferring d = 2, 3, and 4— level states 
through half way distances in closed rings, using this strategy and compare it with 
the original Bose strategy. Thus when r — s = 20, we are using a ring of size N = 40. 

There are a few interesting features. First it is seen that the average fidelity is 
almost the same in both strategies, which implies that we can always, even for d = 2, 
use the second method which is much simpler and does not rely on distance-dependent 
tuning of magnetic fields. Second we note that the strong similarity of the fidelity 
curves in various dimensions. To see the reason of this universality, we note that for 
very small magnetic fields, Bt « 1, where Td(Bt) ~ (d — 1), the average fidelity 
behaves as 
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The curves in figure (1) show the fidelities at the optimal time t op where (F(t opt )) 
becomes a maximum, obtained by numerical searches in a time span t £ (0, 400). Let 
us now suppose that the optimal time is the time where cos(7 rs ) ~ 1. Since f rs is 
independent of d, we can obtain a universal relation for optimal fidelities from the 
above equation by rewriting it as follows. First we note from the above equation that 



which is a constant. To check this assumption and the resulting universality, we draw 
in figure (2) the left hand side of (|7|) (as obtained from numerical searches for the 
optimal time, leading to figure (1) and not by setting cos7 rs « 1) for several values 
of d. The universal behavior is now completely evident. 

Finally, we note that almost perfect state transfer of any d— level state is possible 
when \r — s\ = 1,2,4. This possibility of almost perfect state transfer was first noticed 
in pQ for d = 2 level states. We now see that this is a general and curious feature 
of random swapping for any d— level state. We know that by sequential swapping 
at any two consecutive sites, one can perfectly transfer an unknown state through a 
chain. However this requires a multitude of control operations at all sites of the chain. 
The above result about perfect transfer over distances of 4 sites by random swapping 
(induced by the natural Hamiltonian dynamics), means that one can transfer states 
perfectly over long distances in a chain by a smaller number of control operations, 
namely by 1/4 of the number of sites of the chain. 

Figure (3) shows average fidelities for three different distances, namely |r — s\ = 
4, 7, 14, as functions of the number of levels d. It is seen that the average fidelity 
decreases with d and saturates to a constant value, depending on the distance. Plots 
a and b refer to two different strategies discussed above. 

3 Entanglement distribution 

One of the major problems of quantum information science and technology is the 
distribution of entangled pairs over long distances. For flying qubits, such pairs in 
the form of polarization-entangled photons have been distributed to various long dis- 
tances through optical fibres and free air [201 1211 [22l 1231 [24, 25j. For small distances in 
a quantum computer, for which we supposedly will be dealing with solid state devices 
or ion traps in the future, one needs to distribute entangled pairs through such chains 
of qubits. In [lj a method was proposed for such a task, when the particles have spin 
1/2 or have two levels (qubits). Here we generalize this idea to d— level states or 
qudits. We will see here that the quality of entanglement transfer is better for higher 
values of d. 



F* opt := lim (F opt ) = \f r , 

a >oo 



Rearranging the terms of ([6]) after setting cos 7, 



rs 



1, we find 
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Figure 1: (Color online) The average fidelity for transferring d=2, 3 and 4— level 
states for two different strategies explained in the text. In each case the distance 
between the sender (s) and the receiver (r) is half the length of the chain (N). 
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Figure 2: (Color online) The universal relation for optimal fidelities, explained in 
equation (J7J). 
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Figure 3: (Color online) The average fidelity for three different distances, namely 
|r — s\ = 4, 7, 14, as functions of the number of levels d. Plots a are for the strategy 
of Bose pQ and plots b are for our strategy. 
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Suppose that a maximally entangled state is prepared between sites (not coupled 
to the chain) and site s of the chain. We want to use the natural dynamics of the 
chain to transfer this entanglement through the chain. In particular we want to see 
what will be the entanglement between sites and r at time t. The initial states is 

^ d-l 

\iPme)o, s = -7=53 ( 8 ) 

At time t, the joint state of sites and r can be determined by the extension of the 
map ([3]) to one acting on the chain and the external site 0: 

d-l 

P0,r(t)> = Y,( I ® A h)Wme)(4>Me\(I®AI). (9) 
fj.=0 

Insertion of the Kraus operators A^s from ©in ([9]) gives, 

A)r(i) = i{ |00)(00| + X:^;|00)(H + E/-IM( 00 l 

fj.=l A* = l 

+ E /r M X*M(H + E(i - |/^| 2 )|^0)(/,0| }. 

/i,^=l £i = l 

We use Logarithmic Negativity (LN) as a measure of entanglement of the state 
po r which is defined as, 



LN(p 12 ) = log 2 (\\p%\\), \\0\\ = trV&O, (10) 

where by the superscript T 2 , the partial trace over the second space is implied. Log- 
arithmic negativity is an entanglement monotone which is additive and does not 
increase on the average under all partial transpose preserving operations [27]. For a 
pure maximally entangled state like ([8|), equation (fTU|) yields LN(\iJ)me}) = log 2 d. 
In order to calculate logarithmic negativity we need eigenvalues of p Q r r p Q T r . Straight- 
forward calculations shows that 

d-l d-l 



' {|<K>M<K>| + 

fj.=l M = l 



p^pS = ^{|oo)(oo| + El/r At s lVo)^o| + El^l 2 |o^)(o/.| 



+ E i/^i 2 i/;j 2 i^xh + E(i-i/^i 2 )Vo)^oi 

(J,,V=l M = l 

+ E(i - \a 2 )f? s \vo)(op\ + E(i - \f? s \ 2 )f?;\op)(po\}. 

fj.=l A* = l 

It is readily seen that this operator is the direct sum of two-dimensional matrices of 
the form 

]_( \m\ 2 (i-i/ai 2 )/£ \ 
<p { (i-\m\ 2 )m i-\m\ 2 +\m\ 4 
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1 I f I 

in the subspaces spanned by |^0) and |0/i) (with eigenvalues Ur Jy ) and a diagonal 
matrix spanned by the rest of basis vectors. Putting these together we find the 

spectrum of the matrix \J Pq^ Par as follows, 

i d 

I I Iff. I ( A i\2 



(d-1) 
d-1, 



where the number in front of each eigenvalue denotes its degeneracy. Therefore the 
logarithmic negativity of the final state between sites and r can be computed easily. 
It is found that 

LN( P0r (t)) = iog 2 {i + ± \mfrs\ + \T,\m 2 } 

and because \f{? s \ is independent of p so we can simplify the above formula,LA/'( / oor(i)) = 
log 2 {l + \f rs \ 2 (d — 1)}. This equation shows that with increasing d the logarithmic 
negativity and hence the entanglement increases and indeed approaches its maximum 
value for continuous variable states. We can define the efficiency of entanglement 
distribution as a measure of the the percentage of entanglement that is gained after 
distribution of the maximally entangled state through the chain, so we introduce the 
efficiency as, 

Figure (4) shows the efficiency of entanglement of sites 1 and 30 in a ring of size 
N = 60 as a function of time for three different values of d. It is seen that the optimal 
time for picking up the states is independent of d and the efficiency is increased by 
increasing the dimension d. 



4 Summary 

We have generalized the protocol of [1] for quantum stater state transfer of qubits to 
transfer of d-level states. On the theoretical side, we can consider the results of [1] 
and the present paper as an answer to the question "with what fidelity can a quantum 
state be transferred through a chain if we use random swapping instead of sequential 
swapping?". The latter method is known to achieve unit fidelity but requires local 
control at every site of the chain. We have shown that 1- the fidelity decreases with 
the dimension d, but reaches a saturated value depending on the distance, and 2- 
that when the sender and the receiver are 4 sites apart, nearly perfect transfer is 
possible for any dimension d. A by-product of our study is that we have proposed a 
much simpler method for state transfer, one in which the magnetic field is kept to a 
vanishingly small value, instead of tuning it to a distance-dependent value as in the 
original protocol of pQ. 
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Figure 4: (Color online) Entanglement transition through a chain. The efficiency of 
transmission of maximally entangled pairs through a distance of 30 sites in a chain of 
length 60 for different d— level states. The curves from bottom to top correspond to 
d = 2, d = 3 and d = 4. 

Furthermore the concept of entanglement distribution has been studied for d-level 
states and the interesting result is that the quality of entanglement distribution will 
be improved by increasing the dimension d. 
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